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A GRAPHICAL METHOD OF DETERMINING PRESSURE DISTRIBUTION IN 

TWO-DIMENSIONAL FLOW 

By Robert T. Jones and Dobtb Cohen 



SUMMARY 

By a generalization of the Joukowski method, a pro- 
cedure is developed for effecting localized modifications of 
airfoil shapes and for determining graphically the resultant 
changes in the pressure distribution. The application of 
the procedure to the determination of the pressure distri- 
bution over airfoils of original design is demonstrated. 
Formulas for the lift, the moment, and the aerodynamic 
center are also given. 

INTRODUCTION 

It is possible, by a simple geometric construction, 
to modify any given streamline sbape in sucb a way 
tbat the effect of the modification on the pressure dis- 
tribution, tho lift, and the moment can be readily 
determined. Tho construction is essentially that used 
in deriving the familiar Joukowski airfoils from a circle 
(reference 1) although it may be applied to an airfoil 
shape to introduce modifications of the outline. The 
method is based on the concept of complex numbers 
but its application requires no familiarity with them. 

By two or more successive applications of the con- 
struction to a circle, it is possible to derive shapes of 
such diversity as to permit the approximation of nearly 
any airfoil of current design. For the airfoil derived 
in this way the pressure distribution can be determined 
exactly, so that the only error likely to occur is in the 
reproduction of the exact airfoilshape. Any inaccuracy 
in this approximation is immediately apparent and can 
be made as small as is considered desirable. 

MODIFICATION OP AN AIRFOEL 

The method consists entirely in applications of a 

single construction; desired effects are obtained by the 

proper choice of a parameter k and of the location of 

the axes with respect to the figure to be transformed. 

This basic construction applied to an airfoil surface is 

demonstrated in figure 1. In order to find the point 

z a on the modified airfoil to correspond to a point z 

on the original airfoil, the vector Oz is drawn from the 
u 

origin and the vector 0— is added to it. The vector 

z 

k 2 

0—, which will be termed the "reciprocal" vector, is con- 
structed with its length equal to k* times the reciprocal 
of the length of Oz and at an angle — 0 with thez-axis, 
where is the angle made with this axis by Oz. A 



point of the modified airfoil, then, is located by the 
resultant of the vector Oz and its reciprocal vector. 

The construction by which a figure is modified ap- 
plies equally well to all its streamlines, altering them 
to conform 'to the distortion. The change in spacing 
of the streamlines near the boundary of the figure will 
show directly the effect of the transformation on the 
velocity in that region because the velocity varies 
inversely with the spacing. The factor by which the 
elements of length near a point z are changed may be 
shown to be the ratio of the length of the diagonal from 




Figure i.— Construction to determine a point on tho modified airfoil. 

P/z to z (the vector z—Pfz) to the vector Oz. Then the 
ratio of the velocity at a point z+k*/z to the velocity 
at the corresponding point z of the original figure is 
the inverse of this factor. 

The problem of deteimining the pressure distribu- 
tion is reduced by Bernoulli's relation P=P 0 —pV i f2 
to that of finding the distribution of the velocity. If 
the velocity distribution over the original figure is 
known, the velocity distribution over the modified 

figure may then be found by applying the ratio , t z l 

z 

I z 

at euch point used in the construction. The vectors 
z and z—k 2 /z have already been used in transforming 
these points and are therefore directly measurable. 

Strictly speaking, the method as outlined is applica- 
ble only to potential flows. It is reasonable to suppose, 
however, that the actual velocity distribution over the 
modified airfoil may be obtained with good accuracy 
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from a distribution experimentally determined for the 
original airfoil, provided that: (1) the experimental and 
the theoretical distributions do not differ too greatly 
and (2) the modification itself does not introduce too 
great a change in the distribution. 

The requirements of each problem will suggest the 
proper choice of the axes and of h. A few general 
observations may be helpful in this connection. It is 
evident that the origin must be in the neighborhood of 




Figure 2.— -The effect of varying the parameters ol the construction. 



the section to be modified since the point nearest the 
origin is shifted the greatest distance. The direction 
in which a point ia shifted depends on the inclination of 
its vector.to the axes, and the shift may have a relatively 
small component normal to the airfoil boundary. If 
one axis is roughly parallel to the section of the surface 
to be modified, the point of greatest deviation will 
usually be very near the. intersection of the surface 
with the other axis. It will therefore be useful to note 
that the displacement of a point on either axis is 
k*fd, where d is the distance of the.origin from the point. 
Points on or near the z-axis are shifted outward from 
the origin; points near the.y-axis are moved inward. 
The transformations of figure 2 show the nature of the 
modifications to be obtained by various choices of the 
axes and of k. Sharp modifications occur when the 
airfoil intersects either axis at a distance from the 
origin only slightly greater than k. 



It is expected that the method as outlined will be 
useful in ascertaining the effect on the velocity distri- 
bution of a localized modification of an airfoil. It is 
possible, however, to approximate a more extensive 
modification to any desired degree of. accuracy by 
successive applications of the same transformation, the 
velocity being calculated at each stage by the rule 
already given. 

APPROXIMATION OF A GIVEN AIRFOIL 

The foregoing discussion is concerned with the 

problem of effecting small modifications of existent air- 
foils for which the velocity distribution is known. It 
is sometimes required to predict the theoretical charac- 
teristics of an airfoil not derived in this way. In sucli .. 
a case; it is customary to use the known Sow around a 
circle as a starting point. If the transformation 
discussed here is applied to a circle, the result is an 
oval shape with circular-arc camber, as demonstrated 
by Joukowski. (See reference 1 .) It should be possible 
to proceed to modify this figure (or its special form, a 
Joukowski airfoil) as was done in the preceding section 
to an_experimentally known airfoil. In general, how- 
ever, it would require many steps to reproduce an 
arbitrary airfoil in this way because modifications of 
the shape already approaching an airfoil would, of 
necessity, be small and localized. On the other hand, 
if the procedure were reversed, it would be found that . 
any airfoil may be derived by a single step from a 
figure closely approximating a circle, a figure which 
will hereinafter be called the "distorted circle" of the 
airfoil. The distorted circle may then be considered 
the result of modifying a circle, a result obtained in 
the same way as were the slightly modified airfoils of 
the preceding section. . 

For the first step in this process, the determination 
of the distorted circle corresponding to the desired 
airfoil, reference could be made to Theodorsen and 
Garrick (reference 2) or von Karman" and Burgers 
(reference 3), who give exact formulas for the distorted 
circle in terms of the airfoil coordinates. A graphical 
method based. on the simpler transformation that is 
being used in this paper makes it possible, however, to 
obtain the distorted circle merely by trial. 

The entire procedure is illustrated step by step in 
figured. In order to achieve a considerable simplifica- 
tion of the construction, the airfoil is considered in this 
figure and in the following discussion to have been 
drawn through the midpoints rather than through the- 
tenmnals of the. vectors z+W/z, so that they appear 
here as if to half scale. 

The axes and the parameter k for the first step are 
found from the airfoil dimensions according to the 
relations given in figure 3(a). The z-axis is made to 
pass through the. leading and the trailing edges to 
reduce to a minimum the distortion that will later have- 
to be reproduced. 

The axes and k having been chosen, the intersections 
of the distorted circle with the axes, X N , X T , Y v , and 





Fiqtoe 3.— Steps in the approximation of a given airfoil by the transformation of a circle. 
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i'/,, arc easily calculated (fig. 3(b)). The construction 
then proceeds as follows: 

Id order to serve as a first approximation of the 
distorted .circle, a circle is. drawn to fit the four inter- 
sections with the axes as closely as possible. It should 
be remembered that, because this step is a first approxi- 
mation, great accuracy is not a consideration. When 
the vectors z follow tins circle, the vectors P/z, con- 
structed as previously described, are found to terminate 
in a smaller circle, which is termed for convenience the 
"reciprocal" circle. This reciprocal circle may be easily 
located from its intersections with the axes (fig. 3 (c)). 



Vectors Oz _and 0— are then drawn from the origin 
z 

to the circle and to the reciprocal circle at equal angles 
below and above the z-axis (fig. 3(d)). The midpoint 
z a ' of the diagonal z—tflz connecting the ends of these 
vectors will not, in general, fall exactly on the airfoil. 
A point of the distorted circle z» is located by adjusting 
with dividers the length of Oz b along Oz so that the 
diagonal to its reciprocal is bisected, by the airfoil. 
The point of the airfoil z a at which this bisection occurs 
mil actually be indistinguishable from the intersection 
of the first diagonal with the airfoil; this fact is useful 
in making the correction. Because the corrections to 
k 2 

the vectors Oz and 0— are in approximately the same 

ratio as the vectors, the decrease in k*jz that corres- 
ponds to an increase in the length of Oz may be esti- 
mated and the reciprocal relation be maintained by 
eye. Any desired degree of accuracy may be obtained 
by checking numerically the lengths of the vectors. 

The distorted circle obtainod in this way must now 
bo approximated by applying a transformation of the 
form z+k 2 lz to the circle most nearly approximating it. 
For this purpose the circle already drawn for the pre- 
ceding step will usually be satisfactory, although it 
may havo to bo shifted slightly to obtain a more con- 
venient distribution of the distortion. Depressions 
should occur as nearly as possible on opposite sides of 
the circlo, with distended portions betwoen them. 
(See fig. 3(c).) The ar-axis should then be passed 
tlirough the distended portions and the y-axis, through 
the depressions. Some adjustment will be needed to 
satisfy the condition that the segmonts of the axes cut 
off by the circle and the corresponding distortions be 
approximately in inverse proportion, that is, that thoir 
product at each intersection be constant. The param- 
eter k for the transformation is then tho square root 
of this product. 

Carrying out the transformation of the circle at 
this point servos at once to check the estimated param- 
eters or to suggest an adjustment and to provide the 
construction lines that will be needed to determine tho 
velocity distribution. Similarly, the distorted circle 
thus obtained from tho circle should be reduced to the 
desired airfoil. 

In figure 3(f) the complete double transformation 
and the lines necessary for the computation of the veloc- 



ity distribution arc shown applied to a point of tho 
circle of figure 3(e). 

Occasionally an airfoil will givo riso to a distorted 
circle that cannot be obtained by a single modification 
of a circle. In such a case an additional transforma- 
tion may be applied either to achieve the desired 
distortion of the circle or to modify locally an airfoil 
that can be derived by two transformations. 

VELOCITY DISTRIBUTION OVER AN AIRFOIL DERIVED 
PROM A CIRCLE 

It is apparent that a large variety of .useful airfoil 
shapes can bo obtained by two or more successive appli- 
cations to a circle of the [transformation 'z-{-k 2 jz. 
(Two typical examplos are shown in fig. 4.) The first 
transformation primarily determines the general outline 
of tho airfoil. The second transformation reduces the 
figure to the dimensions of an airfoil and determines 
the nose radius, the thickness, and the camber. It 
may be of interest to note that in the first transforma- 
tion k is small relative to the radius of the circle; in 
the second transformation the ratio of k to the radius 
of tho circle is only slightly less than 1. 

The method of finding the velocity at a point on an 
airfoil derived in this manner is deduced from tho 
following considerations: 

The velocity at a point z on the circle is given by 
the formula 

V,=2V t sin 

where r is the radius of the circle, F 0 is the wind veloc- 
ity, and Q is the angle that the radius to the point z 
makes with the direction of the air stream. The circu- 
lation T is determined by the Kutta condition. Then 

F,=2F 0 (sin e-sin 6 0 ) 

where 9 0 equals the angle between the air stream and 
the radius to z T j., the point of the circle that trans- 
forms to the trailing edge of the airfoil. This expres- 
sion for V t lends itself readily to graphical evaluation 
as a part of the construction; the velocity factor 
(sin 8— sin 0 O ) is proportional to the ordinate of the 
circle measured from a line drawn through z T . B . and 
parallel to the wind velocity V 0 . 

The velocity at the corresponding point of the trans- 
formed figure has been shown in a preceding section to 



be 



¥ 
z— 



V t . Tho "stretching factor" 



k-- 



will have 



to be applied again to the velocities over the distorted 
circle to transform the flow to that over the airfoil. 
A simple procedure is to plot the stretching factor alone 
for the first transformation, as a function of the angular 
position of the transformed point, with the lengtlis 



and 



2 measured from the construction. 

2 . . 



The 



stretching factor for the second transformation could 
then be- applied at convenient points to values of the 
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Figcrz 4.— Airfoils derived by two applications of transformation to circle. Airfoils are shown half scale with respect to circles. 



stretching factor for the first transformation read from 
the plot. The resultant stretehing factor obtained in 
this way can be directly applied to the velocity at 
points of the circle to give the velocity distribution 
over the airfoil at any angle of attack. 

Figure 5 shows a sample pressure distribution ovei 
au airfoil (the NACA 4412 airfoil) approximated by 
two transformations of a circle. Comparison is made 
with the distribution derived by the theory of reference 
2 and with experimental results taken from reference 4. 
The theoretical distributions were computed at an 
angle of attack of 6.4°, which corresponds to a geomet- 
ric angle of attack of 8.5° for the finite-span airfoil of 
reference 4. 

CONCLUSION 

The method of this paper has been found useful in 
determining the effect of small modifications of airfoil 
shapes on the pressure distribution. It is apparent 
that a family of related shapes can be derived in this 
way with greater simplicity than by standard methods 
because the effect of the modification alone can be 
calculated for each shape. It is also possible to foresee 
the maimer in which a shape- must be modified to 
produce a desired change in the pressure distribution. 
The method is, in fact, reversible and by it an airfoil 
may be designed to have a predetermined pressure 




Fiovee 5— Tbeoa'tical and experimental pressure distributions over NACA 4412 
airfoil. ct.=G.4°. 



distribution, provided that a somewhat similar airfoil 
is already known. The modified pressure distributions 
obtained in this way have closely checked with ex- 
periment. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., November 12, 19A0. 



APPENDIX 



FORMULAS FOR LIFT AND MOMENT 



In the ease of an airfoil derived from a circle, the 
method of derivation, interpreted in terms of complex 
numbers, makes it possible to give concise formulas for 
the lift and the moment in terms of the parameters 
used in the construction and to locate by simple geo- 
metric methods the theoretical aerodynamic center. 
The general theory for a transformation, of the form 



Z <L 



(1) 



has been developed by von Mises (reference 5, ch. VII). 
In order to apply the formulas of von Mises, it is neces- 
sary to express the resultant of two successive trans- 
formations in a series of the foregoing form. 

Let z b =z-\rk?l2 express the first transformation, z 
being a point on the circle or one of its streamlines 
and z bl the corresponding point associated with the 
distorted circle. The axes for the second transforma- 
tion are described bj T the complex parameter z u which 
locates the origin, and /3, the angle between the two sets 
of axes. 

Then, if primes denote the vectors drawn to the new 
origin, 

3,'= (*»-*,)«-¥ (2) 

Substitute 



&=2+ 



Then 



(3) 



which is the point on the distorted circle located with 
respect to the axes for the second transformation. 
Apply the second transformation to z b ' ; then 



(4) 



gives the corresponding point associated with the air- 
foil. Substitute for z b ' from equation (3). Equation 
(4) then becomes 



In order to restore the wind velocity T 7 0 to its original 
magnitude and direction, it is necessary to return to 
the original axes. 
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z a '={z a — zi)e~^ 



and the inverse transformation is 



Then 



Za = Z- 



Z 



'^—ZiZ+ki* 



(6) 

m 

(8) 



This expression for the complete transformation can be 
expanded, by carrying out the division of the last term, 
into the series 



3=2+ — _ 1 — ___(. 



(9) 



(10) 



which is in the form of equation (1), where 

02=Wg, t 

The formulas for the lift and the moment as given by 
Glauert (reference 5, pp. 84 and 85) may now be applied. 
The circulation, and consequently the lift, is unchanged 
by the transformation. Thus the lift depends only on 
the radius of the original circle. 



L=pV a T=8icrp T 0 t sin 0 O 



or 4*t|T'V sin 0„ for the half-scale airfoil 



(11) 
The value 

of ilio, the moment about 0 (fig. 6), is given by the 
imaginary part of the expression (reference 5, p. 84) 

r 2- 



) rp^2a 1 TVV'«-2r 2 lV- 



TT 4 A 



3 



(12) 



where the quantities not already defined are as defined 
by the figure. Substitution for a-i from equation (10) 
gives 



TP 



[2 (y+y^e^-^ie""- 



for expression (12), or 

M 0 =27rpV 0 2 [fci 2 sin 2a+k 2 * sin 2(a+0)] 

—mpVoT cos {a— S) 

Since 

r=4?rrT 7 o sin (a— «o) 

Mo=27rpF 0 2 [ti 2 sin 2a+Zr s 2 sin 2(a+(3) 

— 2rm cos (a— 5) sin (a— «o)] 



(13) 

(14) 
(15) 

(16) 
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Figure 6. -Symbols used in development of lift and moment formulas. 



In general, the moment about any point z (see fig. 6) 
is given by 



M^xpVfth 1 sin 2a+k 1 i sin 2(a+/3) 
+2rd cos (a— 4>) sin (a— ao)] 



(17) 



or 



o^^=[fci ! +fe 2 cos 2/S+rd cos (*+<*,)] sin 2a 
+ [i" 2 2 sin 2/3— rd sin (oi+ao)] cos 2a 
+rrfsin (ai-ao) (18) 

The location of the aerodynamic center (or focus) is 
determined by the condition that the moment about 
that point be independent of the angle of attack. In 
order to satisfy this condition, the coefficients of 
cos 2a and of sin 2a in equation (18) must vanish 
simultaneously. 
Thus 

-fa'+h' cos 2/3) = — <Z cos (0+ao) 



_. 



and 



1^ sin 2/3=<2 sin (o>+ao) 
r 

(e3+a 0 ) = tan-> ~ k » « 2 " 



k*+ki* cos 2/3 
rf i =J(IV+2W cos 2/9+*,*) 



(19) 
(20) 

(21) 



Equations (21) are found to have a simple geometric 
representation. If k*/r and kf/r are two sides of a 
trianglo and (180°— 20) is the angle between them, then 
d is the length of the third side and (o3+ao) is the 
supplement of the angle opposite kf/r. This triangle 
may be used to locate the aerodynamic center directly 
on the construction, as shown in figure 7. As in the 
preceding illustrations, this airfoil has been drawn to 
half scale with respect to the circle and the distance d 
has therefore been bisected. The moment about the 
aerodynamic center 




Fiul'be 7.— Construction to find the aerodynamic center. 

M 0 =2xpTVnf sin W>-«o) (22) 

is obtained directly from equation (IS). The moment 
for the half-scale airfoil is one-quarter of this value. 
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